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Abstract 



We have calculated the atomic electric dipole moments (EDMs) induced 
in 199 Hg, 129 Xe, 223 Rn, 225 Ra, and 239 Pu by their respective nuclear Schiff 
moments S. The results are (in units 10 _17 5'(e fm 3 )~ 1 e cm): d( 199 Hg) = —2.8, 
c2( 129 Xe) = 0.38, d( 223 Rn) = 3.3, d( 225 Ra) = -8.5, d( 239 Pu) = -11. We have 
also calculated corrections to the parity- and time-invariance-violating (P, T- 
odd) spin-axis interaction constant in TIF. These results are important for 
the interpretation of atomic and molecular experiments on EDMs in terms of 
fundamental P, T-odd parameters. 



PACS: 32.80.Ys,31.15.Ar,21.10.Ky 
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I. INTRODUCTION 



Recently, a stringent upper limit on the ground-state atomic electric dipole moment 
(EDM) of 199 Hg was obtained [|, 

d( 199 Hg) = -(1.06 ± 0.49 ± 0.40) x 10" 28 e cm . (1) 

(The respective errors are statistical and systematic.) This is the best experimental upper 
limit on an atomic EDM to date. Combined with calculations, this limit can be interpreted in 
terms of limits on fundamental parity- and time-invariance- violating (P, T-odd) parameters. 
These limits tightly constrain competing theories of CP violation. 

Hg has closed electronic subshells, J = 0. The measurement (|IJ) is therefore sensitive to 
P, T- violating mechanisms which originate from the nucleus. The P, T-odd nuclear moment 
that can induce a Hg EDM is the Schiff moment (the nuclear EDM is screened by atomic 
electrons @ and the magnetic quadrupole moment does not contribute due to zero electron 
angular momentum). Note that there are other mechanisms by which the Hg EDM can be 
induced, such as the P, T-odd electron-nucleon interaction (see, e.g. 0), and the electron 
EDM (it contributes due to the hyperfme interaction) [f|||. 

A value for the EDM of Hg induced by the Schiff moment S was estimated in ||, 
c/( 199 Hg) = —4 x 10 -17 S'(e fm 3 ) _1 e cm . This value has been used for the interpretation of 
the measurements of the Hg EDM in terms of P, T-odd nuclear parameters. However, this 
value was obtained indirectly from an atomic calculation of the EDM of Hg induced by 
the P, T-odd electron-nucleon tensor interaction. 

In this work we have performed numerical calculations of the EDMs induced by S in Hg 
and in other diamagnetic atoms of current experimental interest. Our result for Hg, 

d( 199 Hg) = -2.8 x 10- 17 (— ?V)e cm , (2) 

v e fm 7 

differs from the previous estimate (about 40% smaller) and places a more conservative 
constraint on S, and hence the underlying P, T-odd mechanisms that induce it. 

The other atoms we have studied in this work are 129 Xe, 223 Rn, 225 Ra, and 239 Pu, in their 
ground-states. There has been a recent measurement of the EDM induced in 129 Xe || , 

rf( 129 Xe) = (0.7 ± 3.3 ± 0.1) x 10~ 27 e cm . (3) 

A calculation for the Xe EDM induced by S has previously been performed at the Hartree- 
Fock level 0, rf( 129 Xe) = 0.27 x 10 _17 5'(e fm 3 )"^ cm . In this work we have found that 
with core polarization included this value becomes 

rf( 129 Xe) = 0.38 x 10- 17 (— =V)e cm . (4) 

v e fm ' 

This value is 40% larger than the previous calculation and hence places a tighter constraint 
on the P, T-odd parameters extracted from the Xe measurement. 

No experiments or direct calculations have been performed for the ground-state EDMs 
of Rn, Ra, and Pu. In these atoms there is a possibility for an enhanced Schiff moment due 



to static |T0] or even soft |TTJ octupole deformation. 
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Measurements of EDMs of Rn isotopes are in preparation |L2 |. Rn is a heavier analog 
of Xe. As well as the nuclear enhancement in Rn compared to Xe, there is an electronic 
enhancement arising from the increase of the electronic density near the nucleus with higher 
Z . An estimate of the Rn EDM was made in Ref. JlO| by extrapolation of the EDM of Xe 
||, taking into account the enhancement due to increase in Z . 

Also, estimates for Ra [jTUJ] and Pu ]TJ|] have previously been carried out using the same 
extrapolation method as for Rn, although here the results were extrapolated from the esti- 
mate of the EDM of Hg §. 

In Section [TV] (Table |V11] ) the values of previous calculations/estimates for the atomic 
EDMs induced by S in Xe, Rn, Ra, Hg, and Pu and the results of this work are presented. 

The experimental study of the P, T-odd effects induced by the nuclear Schiff moments is 
not restricted to atoms. In the paper JTJJ the P, T-odd spin-axis interaction in TIF molecule 
was measured. From this experiment a limit can be placed on the Schiff moment of the 205 T1 
nucleus. Recently the accuracy of the molecular calculations was significantly increased (see 



15| and references therein) and here we want to use these new molecular calculations to 
obtain a more reliable limit on S^^Tl). This is dealt with in Section ptTJ. 



II. METHOD OF CALCULATION 

The nuclear Schiff moment S produces a P, T-odd electrostatic potential <p. The inter- 
action Hamiltonian 

Hw = J2 h w = -^T, (5) 

i i 

mixes states of opposite parity and induces a static atomic EDM. In previous calculations 
the contact form for the electrostatic potential was used, 

(p(R) = 4ttS ■ V5(R) . (6) 

However, for relativistic electrons the matrix element of this potential diverges. Usually this 
problem is solved by a cutoff of the electron wave functions at the nuclear surface. In the 
work [JT^] it was found that there is a more convenient form for ip suitable for relativistic 
atomic calculations, 

y>(R) = --^p(R) , (7) 

where B = J p(R)R 4 dR, and p(R) is the nuclear density. This expression arises because the 
Schiff moment produces a constant electric field, along the nuclear spin, inside the nucleus 
0- 

The atomic EDM induced in the many-body state N by the P, T-odd interaction Hw 
can be expressed as 

d _ 2 ^ (N\H w \M)(M\D z \N) ^ (g) 
M En — Em 
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where the sum M runs over a complete set of many-body states, En and Em are atomic 
energies, and D z is the atomic electric dipole operator. 

The atoms we have studied in this work are 129 Xe, 223 Rn, 225 Ra, 199 Hg, and 239 Pu. Of 
course, the method of calculation of the electronic states depends on the electron configu- 
ration of the atom. The starting point for the calculations is to self-consistently solve the 
single-particle relativistic Hartree-Fock (HF) equation 

ho\i) = €i\i) (9) 
ho = ca ■ p + (/? - l)c 2 - Za/r + V (10) 

for each electron \i) contributing to the potential V. 

For the noble gases Xe and Rn the calculations are performed in the Hartree-Fock ap- 
proximation in the V N potential. This corresponds to solving (^,0) for the N electrons of 
the atom, that is, in a self-consistent potential V N formed from all N electrons. 

The atoms Ra and Hg can be treated as closed shell systems or as atoms with two 
electrons above closed shells. Correspondingly, we have performed two separate calculations 
for these atoms: one in the V N potential; and a more detailed calculation in the V N ~ 2 
potential (solving (P,|10D for the N — 2 electrons in the atom) with correlations between the 
external electrons and between the external electrons and the core included. The interaction 
of the two external electrons is treated using the configuration interaction (CI) method 
and the correlations of the external electrons with the core accounted for by many-body 
perturbation theory (MBPT). This technique, the CI+MBPT method, was developed in 
the works HJJ. 

Pu corresponds to an open shell atom, with electron configuration 5/ 6 7s 2 . We have 
performed a simple calculation for this atom, in the V N approximation, but with the con- 
tribution of the open /-shell weighted according to its occupancy. (This simple calculation 
in the V N potential, with no correlations accounted for, is justified by comparison of the 
results of the calculations for Hg and Ra in the potentials V N and V N ~ 2 ; see Section |111 B| .) 

As a test of our wave functions we have performed calculations of the ionization potentials 
and the scalar polarizabilities a of the ground-state for each atom. These can then be 
compared with the available experimental data. It is easy to calculate the polarizability 

(ID 



a 



M En — Em 



by replacing the operator H w in (H) by the dipole operator D z 
Below we outline the V N and V N ~ 2 calculations. 



A. V approximation 

In the V N calculation, we can write the atomic EDM induced by the Schiff moment as 

d = 2^n w \d z \n) , (12) 

n 

where the sum runs over the core states \n), d z is the single-particle dipole operator, and 
\Snw) denotes the correction to the state \n) due to the P, T-odd Hamiltonian hw The 
correction \5riw) can be expressed as 
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x ^ (a\h w \n) 

\tinw) = Zs Z ( 13 ) 

where |a) corresponds to an excited state. It is found by solving the equation 

(h - e n )\5n w ) = -h w \n) . (14) 

(Of course, calculating the correction to \n) due to the electric dipole (El) field and taking 
the matrix element of the weak Hamiltonian is equivalent.) 

To include polarization of the core due to the fields hw and d z , we need only include the 
polarization due to one field, e.g. by replacing hw in ([L4D by h w = h w + SVw- (The correc- 
tion 5V\y accounts for the change in the self-consistent potential V N due to the modification 
of the wave functions. For more on the time-dependent Hartree-Fock (TDHF) method, or 
random-phase approximation (RPA) with exchange, see, e.g., [PH .) This is because all the 
states |n) belong to the core (that is, they are included in the Hartree-Fock potential), so 
Y^ n (Sn w \d z \n) = J2n($nw\d z + $V d \n). 

The scalar polarizabilities are easily calculated by replacing the correction \Snw) due to 
the P, T-odd field by the correction \8rid) due the El field. 



B. V approximation 

Hg and Ra can be considered as atoms with two valence electrons above closed cores, 
[ls...5<i 10 ] and [ls...5d 10 6s 2 6p 6 ], respectively. The calculations performed in the V N ~ 2 ap- 
proximation correspond to the CI+MBPT method developed in the works p6| , |T7|| . An effec- 



tive Hamiltonian is constructed for the valence electrons which is formed using MBPT for 
the interaction of the valence electrons with the core. In this way the correlations between 
the valence electrons 6s 2 (for Hg) and 7s 2 (Ra) are treated using the CI method, while the 
valence-core correlations are treated using MBPT. (For more on this procedure, please refer 



to 16,1 



Note that calculations in the CI+MBPT method can be performed in a different poten- 
tial, even in V . The final accuracies for V N ~ 2 and V N potentials are comparable; however, 
calculations in V N ~ 2 for Hg and Ra are somewhat simpler — there are no subtraction 



diagrams, see, e.g., Ref. |T6 



For Hg the CI states were formed from the one-electron basis sets Is — 12s, 2p — 13p, 
3d — 14c?, and 4/ — 15/. The core Is — 5d and the states 6s, 6p, 6d, 5f were HF ones, 
while the rest were virtual orbitals. The virtual orbitals were formed using the recurrent 
procedure, similar to that used in [Hi] . Full CI was made for the two valence electrons. Two 



sets of virtual orbitals were constructed to test the saturation of the CI space. The results 
appeared to be very close, so we concluded that saturation was reached. 

MBPT calculations generally require more basic functions for high energies and more 
partial waves, so the number of virtual orbitals was increased to construct the MBPT basis 
set: it included Is - 21s, 2p - 22p, 3d - 23d, Af - 21 /, 5g - lSg, and Qh - I9h. For Ra the 
CI basis set was somewhat larger and included orbitals Is — 15s, 2p — 16p, 3d — 17d, and 
Af — 17 f. The extended basis set for MBPT was similar to that used for Hg. 
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III. RESULTS 



A. Xenon and radon 

The binding energies of the ground state of Xe and Rn performed in the HF approxima- 
tion in the V potential are presented in Table |I[ It is seen that the Hartree-Fock calculation 
serves as a good approximation for the noble gases, with the calculated and experimental 
ionization potentials in disagreement only at the level of 1%. 

The ground-state scalar polarizabilities a are listed in Table Core polarization in- 
creases a only very slightly (about 1%) from the Hartree-Fock value. The value for Xe is in 
perfect agreement with the experimental value. 

Our results for the electric dipole moments induced in Xe and Rn due to the nuclear 
Schiff moment are presented in Table |ITT| . The effect of core polarization increases the value 
of the EDMs in Xe and Rn by about 30%. Due to the increase in Z, the EDM induced in 
Rn is ten times larger that that induced in its lighter analog Xe. 

A calculation for the Xe EDM has previously been performed at the HF level, with the 
result d( 129 Xe) = 0.27 x 10~ 17 S(e fm 3 )"^ cm §. Our HF result is in agreement with this 
value. 

No direct atomic calculation has previously been performed for the atomic EDM of Rn, 
or for any of the remaining atoms in this work. The estimate for Rn, based on extrapolation 
of the HF result from Xe, is ( 223 Rn) = 2.0 x 10- 17 S(e fm 3 )"^ cm 0. This value is not 
significantly different from our HF value. 



B. Mercury and radium 

The ionization potentials of Hg and Ra performed in the V approximation are listed 
in Table p] alongside the results for the noble gases. The deviation from experiment for Hg 
and Ra is ~ —14%, indicating the importance of correlations. 

The ground-state polarizabilities for Hg and Ra performed in the V N approximation 
are listed in Table ||. We see here that core polarization increases a, from the HF value, 
by about 10% for Hg. The effect of core polarization on Ra is more significant, with an 
increase of 45%. For both Hg and Ra, essentially the entire HF and TDHF results for a 
arise from the 6s 2 and 7s 2 contributions, respectively. We see that the result for Hg is 
strongly overestimated (about 30% higher than the experimental value). This indicates that 
the calculations produce wave functions that are very poor at large distances from the 
nucleus for systems with two electrons above closed shells. 

Our results for the EDMs induced in Hg and Ra due to the nuclear Schiff moment, in 
the V N approximation, are presented in Table [TTT| . The effects of core polarization in Hg 
and Ra are dramatic, with the EDM increased by a factor of 2.5 and 4.5, respectively. The 
dominating contributions to the EDMs of Hg and Ra come from 6s 2 and 7s 2 , respectively. 
These contributions are larger than, and of opposite sign to, the contributions from the core. 
The instability of the results motivates us to study these atoms more carefully, taking into 
account valence-valence and valence-core correlations. 

The V N ~ 2 approximation is more appropriate for calculations of two-electron atoms. 
The results of the calculations for the binding energies of relevant states of Hg and Ra are 
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presented in Table |V]. The calculated value for the removal energy of both s-electrons for 
Hg and Ra deviates from experiment by about —10% at the CI level, and is then improved 
to +2% with MBPT included. This is a significant improvement compared to the result. 
The accuracy of the calculations of the energies of other states is not as good. These states 
are not directly relevant to our calculations of the EDMs, however we have presented them 
in the table as an indication of the accuracy of wave functions produced at various stages 
of the CI+MBPT method. At the CI stage of the calculations the energies of the ss and sp 
states are underestimated by 10 — 20%. Note that for sd states the energies are overestimated 
by about the same amount. With MBPT corrections included, the accuracy of the energies 
improves to 2 — 9%. 

Results of the calculations of the ground-state scalar polarizabilities a of Hg and Ra 
in the V N ~ 2 approximation are listed in Table |V[ We see that for Hg at the CI level of 
calculation a is strongly overestimated (by about 60%). However, both the MBPT and 
TDHF corrections reduce the value so that the final result is in excellent agreement with 



experiment (compare 32.99 a.u. with 34 a.u. |22]|). We see that for Ra the MBPT and 
TDHF corrections bring about a cancellation of the CI value similar to what we see for Hg. 

The results of the calculations of the EDM of Hg and Ra in the V N ~ 2 approximation are 
presented in Table (VJ. We see here that for Hg the MBPT and TDHF corrections increase 
the CI approximation by a factor of two. The corrections for Ra are huge, the final value 
3.6 times that of the CI value. 

Even though these corrections are so large, it is certainly interesting that the result for 
<i(Hg) in the V N ~ 2 approximation (with correlations) is in agreement (~ 10%) with the 
simple calculation in the V N approximation; compare the results of Hg and Ra in Table 
|VI| with those in Table HI, This coincidence of the results is an argument in favour of 



the stability of the final values. Also it indicates that the simple V N approximation gives 
accurate results for EDM calculations of closed-shell atoms. From a consideration of the 
results, we take as our final value for the EDM in Hg induced by S, 

d( 199 Hg) = -2.8x 10- 17 (— ^)ecm . (15) 



-e fm 3 - 

We take as our final value for the EDM induced in Ra 

d( 225 Ra) = -8.5 x 10 17 f — S -^]e cm . (16) 

v e fm 7 

We can compare these new values for Hg and Ra with estimates calculated in Refs. 
and fllQfl , respectively. The value for Hg, <i( 199 Hg) = —4 x 10~ 17 S(e fm 3 ) _1 e cm, was obtained 



indirectly from an atomic calculation |7j of the Hg EDM induced by the P, T-odd electron- 
nucleon tensor interaction. The radium estimate, <i( 225 Ra) = —7 x 10 _17 S'(e fm 3 ) _1 e cm, 
was found by extrapolation of the Hg EDM estimate, taking account of the enhancement 
due to higher atomic number Z. 

C. Plutonium 

We performed a simple calculation for Pu. We basically used the V N approximation. 
The unfilled /-shell was accounted for by weighting the corresponding angular coefficients 
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according to the occupation number of the shell. In calculations of the atomic EDM, we can 
expect that Pu behaves in a similar way to Ra and Hg (there is a closed s-shell, 7s 2 , and 
an open /-shell). The /-shell does not contribute to the EDM at the HF level (the /-shell 
does not penetrate the nucleus due to the centrifugal barrier, and so does not contribute to 
the Schiff matrix element). However, it can contribute due to core polarization. Because 
the V N calculations for Ra and Hg turned out to be sufficient, we expect that the same is 
true for Pu. 

Our result for the ionization potential of Pu (see Table |) deviates by -15% from exper- 
iment. The polarizability a is listed in Table [TT| As we saw for Hg and Ra, essentially the 
entire value is due to the 7s 2 contribution. The results for the EDM of Pu are presented in 
Table |T|. These values are essentially due to the contribution of the 7s-electrons, which are 
of opposite sign to the contributions of the core and /-electrons. With core polarization, the 
contribution of the /-electrons amounts only to 10% of the final value. The effect of core 
polarization increases the HF value 4 times. This is similar to what we see in Ra. 

An estimate for the Pu EDM was performed in Ref. |13| by extrapolation from the 
estimate of the Hg EDM, c?( 239 Pu) = —10 x 10~ 17 S(e fm 3 ) _1 e cm. Our calculation is in 
agreement with this estimate. 



D. TIF molecule 

All molecular calculations deal with the following matrix element: 

X = — ^(*o|[V-n,5(R)]|* ), (17) 

where \l/o is the ground state wave function and n is the unit vector along the molecular 
axis. It is clear that this matrix element is related to the contact form of the Schiff moment 
potential given by Eq. @: 

-4vr(^ |S • [V,5(R)]|* > = 6XS ■ n. (18) 

It is not difficult to find the correction coefficient k\ which accounts for the difference between 
the contact form of the Schiff moment interaction (|(J) and the more accurate expression (0). 



The latest molecular calculation [151 accounted for the correlations between valence elec- 
trons, but neglected core-valence correlations. We saw above that the latter appears to 
be very important in atomic calculations. The most important core-valence correlations in 
molecules with one heavy atom are of the same nature as in atoms and can be accounted 
for by an atomic calculation for the heavy atom | 23fl . The valence space in the calculation 



TBI included 5s, bp, 5d, 6s, and 6p electrons, leaving a rather small core [Is 2 . . . 4/ 14 ]. For 
such a compact and rigid core the dominant correlation correction is from the RPA for the 
Schiff potential. We account for this correction by the coefficient k 2 . 

The final expression, which includes the two corrections discussed above, has the form: 

(V \H w \y ) = 6k 1 k 2 XS-n. (19) 

The dominant contribution to the molecular matrix element ( fl9|) comes from atomic matrix 
elements for Tl of the form 



S 



1 2 

(ns\H w \mp z ) = ~{nsx/2\H w \mpx/ 2 ) + - (ns 1/2 \H w \mp 3/2 ) . (20) 

This allows us to approximate the correction coefficients: ki ~ \ki(si/ 2 , px/ 2 ) + ^ki(sy 2 , P3/2) ■ 
Within this method we calculated the coefficients ki for the Tl atom and found: 

fci = 0.89; k 2 = 1.10. (21) 

We see that both corrections are relatively small and, therefore, our simple model calculation 
is justified. Note that the core here is much smaller than in our previous atomic calculations. 
Because of that the RPA correction here is only 10%. Other core- valence correlations tend to 
decrease more rapidly with the energy of the core orbitals, and therefore are less important 
here. Using the result X = 7635 a.u. from ]15| and Eqs. fll9f) and (|2l| ) we get as our final 



answer: 



{q \H w \q ) = 4.47 x 10 4 (S • n) a.u. . (22) 



IV. CONCLUSION 

The most accurate measurement of the P, T-odd spin-axis interaction in the TIF molecule 



14 1 combined with Eq. (p2|) gives the limit: 

S( 205 T1) = (69 ± 111) x 10~ 12 e fm 3 . (23) 



In Table |VII| we present our final results for the atomic EDMs induced by the Schiff 
moments S alongside the values of previous calculations/estimates. 
The limit on the Schiff moment of 129 Xe (Eq. ([3D, Table [VI I| ) is 

S( 129 Xe) = (184 ± 868 ± 26) x 10~ 12 e fm 3 . (24) 

Our result for Hg places new constraints on the fundamental P, T- violating parameters. 
Comparing (|IJ) and (O) we obtain for the limit on the nuclear Schiff moment: 

S( 199 Hg) = (3.8 ± 1.8 ± 1.4) x 10~ 12 e fm 3 . (25) 

Comparison of the calculated and measured values of the ionization potentials and polariz- 
abilities as well as a comparison of the calculations in two completely different approxima- 
tions (in V N and U^ 2 ) indicate that the error of the atomic calculations probably doesn't 
exceed 20%. 

In general, the Schiff moment can be induced from a number of P, T- violating mech- 
anisms: due to a permanent EDM of an unpaired nucleon or due to the P, T-violating 
nucleon-nucleon interaction. Mercury has an unpaired neutron, and so in the shell model 
its Schiff moment can arise due to the EDM of the neutron. The contribution to the Schiff 
moment of Hg from a proton EDM can be estimated by comparing the experimental value 
of the magnetic moment of Hg with that estimated by the nuclear shell model (this allows us 
to estimate the contribution of proton configurations with unpaired spin); see Ref. |J. The 
dominant mechanism for the production of a Schiff moment is the P, T-odd nucleon-nucleon 



9 



interaction rj p4J . Due to this mechanism, the P, T-odd field of the unpaired neutron excites 
core protons which contribute to the Schiff moment ||. The magnitude of this interaction 
is characterized by the dimensionless constant r] np . In Table |V111| we present the limits on 
these parameters extracted from Hg. These are compared in the table with the best limits 
from other experiments. Also presented are limits on P, T-violating parameters at the more 
fundamental level, the P, T-odd pion-nucleon coupling constant g® NN , the P, T-odd QCD 
phase 9, and the chromoelectric dipole moments (CEDMs) and EDMs of quarks. We see 
here that the limits extracted from the Hg measurement are stronger than those extracted 
from direct neutron EDM experiments. 
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TABLES 



TABLE I. Ionization potentials of Xe, Rn, Hg, Ra, and Pu. The calculated values are obtained 
from the Hartree-Fock approximation in the V N potential. The experimental values are presented 

in the last column. Units: cm -1 . 

Atom HF Exp. a 

Xe 96525 97834.4 

Rn 84285 86692.5 

Hg 71996 84184.1 

Ra 36485 42577.35 

Pu 41463 48890(200) b 



a From Ref. [2C] unless otherwise stated. 
b Ref. pi. 



TABLE II. Scalar polarizabilities of Xe, Rn, Hg, Ra, and Pu obtained in HF and TDHF 
approximations in V . Experimental values are presented in the last column, (a.u.) 



Atom 



HF 



TDHF 



Exp. 



Xe 
Rn 
Hg 
Ra 
Pu 



26.87 
34.42 
40.91 
204.2 
147.3 



26.97 
35.00 
44.92 
297.0 
201.3 



27.06 
34 



a Ref. 



TABLE III. Electric dipole moments d induced in Xe, Rn, Hg, Ra, and Pu by the nuclear 
Schiff moment S. We present results obtained using the HF and TDHF approximations in the V N 
potential. Units: 10~ 17 (5/(e fm 3 ))e cm . 



Atom HF TDHF 

Xe 0.289 0.378 

Rn 2.47 3.33 

Hg -1.19 -2.97 

Ra -1.85 -8.23 

Pu -2.66 -10.9 
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TABLE IV. Binding energies of low states of Hg and Ra calculated in the V approximation. 
The removal energy for both 6s electrons for Hg and both 7s electrons for Ra is presented in the 
first row of each respective atom. Energies of excited states are presented relative to the ground 

state. Units: cm' 1 . 

Atom 



Ra 



State 




CI 


+MBPT 


Exp. a 


6s 2 




207659 


240912 


235464 


QsQp 


3 P 


29336 


40012 


37645.080 




3 Pi 


31009 


41753 


39412.300 




3 P 2 


34794 


46776 


44042.977 




x Pi 


46142 


55395 


54068.781 


6s7s 


3 S 1 


51460 


65044 


62350.456 




1 Sq 


54187 


67090 


63928.243 


7s 2 


<Sb 


115318 


127248 


124419.66 


7s6d 


3 D, 


15910 


14012 


13715.85 




3 D 2 


16067 


14465 


13993.97 




3 D 3 


16625 


15921 


14707.35 


7s7p 


3 Po 


10424 


14268 


13078.44 




3 Pi 


11289 


15159 


13999.38 




3 P2 


13535 


17937 


16688.54 




'Pi 


18835 


21663 


20715.71 



^ef. 



TABLE V. The scalar polarizabilities a for Hg and Ra calculated in the potential V N ~ 2 . The 
valence and core contributions are separated into different rows and their sum is presented. The 
HF result, with CI included into the valence orbital contribution, is presented in the first column 
of results. In the next columns the Brueckner and TDHF contributions are added, fa.u.) 



Atom 



Contribution 



CI 



+ MBPT 



+ TDHF 



Hg 



Ra 



6s 2 

core 

sum 

7s 2 

core 

sum 



46.54 
8.05 

54.59 
321.4 

13.6 
335.0 



36.06 
8.05 

44.11 
260.7 

13.6 
274.3 



25.69 
7.30 

32.99 
218.2 

11.8 
229.9 
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TABLE VI. The EDMs induced in Hg and Ra calculated in the V N ~ 2 potential. The valence 
and core contributions are separated into different rows and their sum is presented. The HF result, 
with CI included into the valence orbital contribution, is presented in the first column of results. In 
the next columns the Brueckner and TDHF contributions are added. Units: l(T 17 (S/(e fm 3 ))e cm . 

Atom Contribution CI + MBPT + TDHF 

Hg 6? T90 -2.77 -3.04 

core 0.63 0.63 0.34 

sum -1.26 -2.14 -2.70 

Ra 7s 2 -4.03 -6.08 -10.10 

core 1.62 1.62 1.40 

sum -2.41 -4.46 -8.70 



TABLE VII. Final results for the atomic EDMs induced by the nuclear Schiff mo- 
ment. These are compared with the simple calculations/estimates of previous works. Units: 
10- 17 {S/{e fm 3 ))e cm . 



Atom Other works This work 
Xe 0.27 a 0.38 
Rn 2.0 b 3.3 
Hg -4.0 c -2.8 
Ra -7.0 b -8.5 
Pu -10 d -11 



a Ref. ||. 
b Ref. gO]. 
c Ref. @. 
d Ref. 
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TABLE VIII. Limits on P, T- violating parameters in the hadronic sector extracted from 199 Hg 
(Eq. (p5|)) compared with the best limits from other experiments. Errors are experimental. Some 
relevant theoretical works are presented in the last column. 



P, T-violating term Value 



System Exp. Theory 



neutron EDM d n 



proton EDM cL 



rinpi-^PPn-fsn 



9ttNN 



QCD phase 9 



CEDMs d and 
EDMs d of quarks 



(17 ± 8 ± 6) x l(T 2b e cm 

(1.9 ± 5.4) x 10~ 26 e cm 

(2.6 ± 4.0 ± 1.6) x 1(T 26 e cm 

(1.7 ±0.8 ±0.6) x 10~ 24 e cm 
(17 ± 28) x 10~ 24 e cm 

f? np = (2.7 ± 1.3 ± 1.0) x 10- 4 



(3.0 ± 1.4 ± 1.1) x 10" 



-12 



(1.1 ±0.5 ±0.4) x 10" 10 

(1.6 ±4.5) x 10" 10 

(2.2 ± 3.3 ± 1.3) x 10" 10 

e(d d - d u ) = (1.5 ± 0.7 ± 0.6) x lO" 26 e cm 
e{d d + 0.5d u ) + l.3d d - 0.3d u 



(3.5 ±9. 



10 



-26 



(4.7 ± 7.3 ± 2.9) x 10" 2b e cm 



e cm 

26 



199 Hg 

neutron 
neutron 

199 Hg 
TIF 



199 



Bp 



199 Hg 

199 Hg 

neutron 
neutron 



199 



Hg 



neutron 
neutron 
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